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Abstract 

The specific topology of the line centered square lattice (known also as the Lieb lattice) induces 
remarkable spectral properties as the macroscopically degenerated zero energy flat band, the Dirac 
cone in the low energy spectrum, and the peculiar Hofstadter-type spectrum in magnetic field. We 
study here the properties of the finite Lieb lattice with periodic and vanishing boundary conditions. 
We find out the behavior of the flat band induced by disorder and external magnetic and electric 
fields. We show that in the confined Lieb plaquette threaded by a perpendicular magnetic flux 
there are edge states with nontrivial behavior. The specific class of twisted edge states, which 
have alternating chirality, are sensitive to disorder and do not support IQHE, but contribute to 
the longitudinal resistance. The symmetry of the transmittance matrix in the energy range where 
these states are located is revealed. The diamagnetic moments of the bulk and edge states in the 
Dirac-Landau domain, and also of the flat states in crossed magnetic and electric fields are shown. 



PACS numbers: 73.22.-f, 73.23.-b, 71.70.Di, 71.10.Fd 



I. INTRODUCTION 



The interest in the line centered square lattice, known as the 2D Lieb lattice, comes from 
the specific properties induced by its topology. The lattice is characterized by a unit cell 
containing three atoms, and a one-particle energy spectrum showing a three band structure 
with electron-hole symmetry, one of the branches being flat and macroscopically degenerate. 
For the infinite lattice, the three energy bands touch each other at the middle of the spectrum 
(taken as the zero energy ), and the low energy spectrum exhibits a Dirac cone located at 
the point F = (tt, it) in the Brillouin zone. Except for the presence of the flat band, the 
Lieb lattice shows similarities with the honeycomb lattice in what concerns both spectral 
and transport properties. For instance, besides the presence of the Dirac cone, the energy 
spectrum in the presence of the magnetic field shows also a double Hofstadter picture, with 
the typical \/B dependence of the relativistic bands on the magnetic field B pp. The Hall 
resistance of the two systems in the quantum regime behaves alike, but the step between 
consecutive plateaus equals h/e 2 in the Lieb case (instead of h/2e 2 for graphene) because of 
the presence of a single Dirac cone per BZ. An all-angle Klein transmission is proved by the 
relativistic electrons in the Lieb lattice [2]-Ej. 

There are more lattices that support flat bands, however it is specific to the Lieb lattice 
that the band is robust against the magnetic field, while other lattices develop dispersion 
at any B ^ 0. The intrinsic spin-orbit coupling does not affect the flat band either, but 
opens a gap at the touching point T, the Lieb system becoming in this way a quantum spin 
Hall phase [3j [5] . Topological phase transitions driven by different parameters are studied 
in [HI [7]. The zero-energy flat bands became a topic of intense study also for other reasons: 
they may allow for the non-Abelian FQHE [BHTO] or for ferromagnetic order and surface 
superconductivity [TT14T3] . 

In this paper we address the properties of the finite (mesoscopic) Lieb lattice with em- 
phasis on some features of the flat band and of the edge states which are specific to this 
lattice. We adopt the spinless tight-binding approach, and the spectral properties are exam- 
ined under both periodic and vanishing boundary conditions applied to the system described 
in Fig.l. In section II we find that the zero energy flat band exists independently of the 
boundary conditions. It turns out, however, that in the periodic case the band is built up 
only from B- and C- orbitals, while in the other case the A-type orbitals are also involved 
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(see Fig.l). We prove this analytically by calculating the eigenfunction in both situations. In 
this way we also find out that for confined systems (i.e., with vanishing boundary conditions) 
the degeneracy of the flat band equals N ce u + 1 {N ceU is the number of cells of the meso- 
scopic plaquette). We find in section IIIA that for a confined plaquette two levels separate 
from the bunch when a perpendicular magnetic field is applied, such that the degeneracy is 
reduced by 2. This is proved in a perturbative manner for the general case, however it can 
be observed more easily by the use of the toy model consisting of two cells only. 

Next, we study how the flat band degeneracy is lifted by disorder and by an external 
electric field applied in-plane. An exotic result is that the extended states of the disordered 
flat band in the presence of a magnetic field behave according to the orthogonal Wigner- 
Dyson distribution although the unitary distribution is expected. When an electric field 
is applied, the flat band splits in a Stark- Wannier ladder whose structure is analyzed by 
calculating the diamagnetic moments of the states in crossed electric and magnetic fields. 

In section IIIB we study the edge states which fill the gaps of the double Hofstadter 
butterfly when the magnetic field is applied on the confined Lieb plaquette. We identify 
three types of such states. The conventional edge states located between the Bloch-Landau 
bands and also between Dirac-Landau bands (i.e., the relativistic range of the spectrum) 
differ, as expected, in their chirality. Additionally, we detect twisted edge states situated in 
the magnetic gap which protect the zero-energy band, coming in bunches and characterized 
by an oscillating chirality as function of the magnetic field. The twisted edge states show 
remarkable properties: surprisingly, they are not robust to disorder, as the other types of 
edge states are, and does not carry transverse current (i.e., the QHE vanishes in the energy 
range covered by these states). The last property comes from a specific symmetry of the 
transmittance matrix which is discussed in section IV. 

Finally, one has to note that the line centered square lattices are found in nature as 
Cu — 2 [2] planes in cuprate superconductors, and can be engineered as an optical lattice 

[SIS]. 
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II. THE TIGHT-BINDING MODEL FOR THE LIEB LATTICE : PERIODIC VER- 
SUS VANISHING BOUNDARY CONDITIONS 



Our aim is to point out specific aspects of the confined Lieb plaquette from the point 
of view of spectral and transport properties. In order to allow for a comparison we shortly 
describe also the case of the infinite system, with and without magnetic field, although the 
eigenvalue problem is already known from the literature. We remind that the continuous 
model for the infinite Lieb system in perpendicular magnetic field [3J , shows the yB depen- 
dence on the magnetic field of the eigenenergies in the relativistic range. The information 
obtained in the long-wave approximation of the Schrodinger equation, concerning the de- 
pendence on B of the Bloch-Landau or Dirac-Landau bands are recaptured in the spectrum 
of the discrete tight-binding model (Fig. 5a) together with effects coming from the periodic 
lattice and finite edges. 

In this section, starting from the tight-binding Hamiltonian, we built up the eigenfunction 
of the periodic and finite Lieb plaquette, and prove the degeneracy and structure of the zero- 
energy flat band. The crossover from the simple Hofstadter spectrum of the simple square 
lattice to the Lieb spectrum characterized by a double butterfly, magnetic gap and a flat 
band is shown in Fig.3. 
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FIG. 1: (Color online) The Lieb lattice: the unit cell contains three atoms A,B and C; indices 
(n,m) identify the cell; t x ,t y are the hopping integrals along the directions Ox and Oy, respectively; 
a is the lattice constant. 



The Lieb lattice is a 2D square lattice with centered lines as shown in Fig.l. It is 
characterized by three atoms (A, B, C) per unit cell, the connectivity of the atom A being 
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equal to four, while the connectivity of atoms B and C equals two. 

Introducing creation {a^ m , b\ m , c4 m } and annihilation {a nm , b nm , c nm } operators of the 
localized states \A nm >, \B nm >, \C nm > (where (nm) stands for the cell index and the 
letters A, B, C identify the type of atom), the spinless tight-binding Hamiltonian of the 
Lieb lattice in perpendicular magnetic field reads: 



H — X/nm E a a\ im a nm + E b b nm b nm + E c c\ 



cC nm C nm 
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(1) 



where is the flux through the unit cell of the Lieb lattice measured in quantum flux units; 

— * 

we mention that the vector potential has been chosen as A = (—By, 0, 0). 

The presence of a spectral flat band can be noticed already in the simplest case of the 
periodic boundary conditions and vanishing magnetic flux. Assuming that the lattice is 
composed of N* ell = N cells along Ox and N^ ell = M cells along Oy, the Fourier transform 
S n m c nm e l ( kxn+kym ^ (and similarly for all the other operators) yields the 



^k C kx,k y y/NM 

^-representation of the Hamiltonian described by a 3 x 3 matrix: 



H = Y j ( a i bl cl) 

^ \ k k k J 



I E a A*(k x ) A*(k y )\ 



A(k x ) 
A(k y ) 



E b 






E r 



I 



( n \ 



\tk) 



(2) 



where k x = 2irp/N (p = 1,..,N), k y = 2irq/M (q = 1,..,M), and the notations A(k x ) = 
t x (l + e ikx ) , A(k y ) = t y (l + e iky ) has been used. With the choice E a = E b = E c = 0, one 
obtains the following eigenvalues: 



n±(k)= ±vWT^ 
n (k)= o, 



±2 



y/tlco3^k x /2)+tlco3^ky/2), 



(3) 



where Q± are the energies of the upper and lower band, respectively, and fl is the nondis- 
persive (flat) band of the Lieb lattice. The most interesting point in the BZ is the point 
T = (vr, 7r), where in the case of the infinite lattice the three branches are touching each 
other. The expansion of the functions A(k x ) and A(k y ) about this point gives rise to a Dirac 
cone (massless) spectrum: 

«± = ^m+tm- (4) 



On the other hand, the expansion of the same functions about R = (0, 0) shows a parabolic 
dependence: 
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(5) 



*2m x 2m y 

where m x ,m y are effective masses along the two directions. 

Other relevant points in the BZ are M = (vr,0) and (0, ir), which prove to be saddle 
points in the spectrum as it can be noticed also in Fig.2. Above and below the corresponding 
energy E = ±2i (where we considered t x = t y = t) the effective mass exhibits opposite signs 
inducing the change of sign of the Hall effect which is visible in Fig. 15. 

For comparison's sake, we remind that the energy spectrum of the honeycomb lattice 
contains two cones per BZ, and that the saddle point occurs at the energy E = ±t. The 
tight-binding spectrum of the graphene extends over the interval [-3t,3t], while for the Lieb 
lattice the interval is [— 2y2t, 2y2t\. 
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FIG. 2: (Color online) The energy spectrum of the infinite Lieb lattice, (left) The case E a = 
Ef, = E c = when the three bands (two dispersive and one flat) get in contact at k = (ir, ir). 
At low energy the dispersion is linear giving rise to Dirac cones, (right) The staggered case 
E a = 0, = E c = 1 when the spectrum is gapped and the rounding of cones is obvious. 



In order to get supplementary information about the origin of the flat band let us consider 
the staggered case E a — 0, E b — E c — E . Then, a gap is expected in the energy spectrum, 
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and, indeed, the eigenvalues are now [I]: 



Ot(fc) = \ [E ± ^ + 4(| + |A| 2 ) ] , fi (^)=^o. (6) 

The new spectrum is shown in Fig.2b, where one notices the persistence of the flat band, 
which is however shifted to E = E . Since the energy E = E corresponds to the atomic 
level of the orbitals B and C, the result argues that the flat band states are created only 
by this type of orbitals. The gap induced by the staggered arrangement is accompanied by 
the rounding of the cones, that indicates a non-zero effective mass in the low-energy range 
of the two spectral branches, as it can be observed in Fig.2b. 

In what follows we shall calculate the eigenfunctions of the finite Lieb lattice, imposing 
first periodic conditions, and then the vanishing boundary conditions proper to the confined 
plaquette. Let ty? be the eigenfunctions of the Lieb lattice with periodic boundary conditions 
built up as the linear combination: 

*s = «s4i° > +w%\ Q > +^4i° >> ( ? ) 

where the coefficients (3^, 7^ satisfy the equations: 

E a a % + A*(k)/? s +A*{k v )>y Ji = Ea i 
A(k x ) ai +E b p i = Ep % 

A(k x )a % +E, n% = E 1% . (8) 
Then, the functions corresponding to the eigenvalues Qq and Q± in Eq.(3) read: 

nk) = -^=(A*(y&t - A*(M4)|0 > (9) 



=-(±4+ . A(M bU- f A{ky) 4)io >. (io) 

V ; 2 V * v/|A| 2 | + |A| 2 k v/|A| 2 + |A| 2 kJl K J 

In the case of periodic conditions applied to the finite plaquette there are some subtleties 
concerning the band degeneracy which become unimportant in the limit of infinite system. 
It is obvious from Eqs.(9-10) that the three bands come into contact at k = (71", tt), however 
this value of k is allowed only if both N and M are even. In this case the flat band at E — 
is (Nceii + 2)- fold degenerate, otherwise all the three bands are A r ceH -fold degenerate (where 
the number of cells N ce u = NM). 
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The expression of ^f°(k) in Eq.(9) indicates again that the flat band of the periodic lattice 
is composed only from orbitals of the type B and C. On the other hand, we shall see below 
that in the case of vanishing boundary conditions the zero-energy eigenfunction may sit also 
on the A— type sites, and that the degeneracy of the flat band becomes N ce u + 1. 

The periodic boundary conditions can be used in the presence of a uniform perpendicular 
magnetic field for rational values of the magnetic flux <p = p/q resulting in a spectrum 
composed of two Hofstadter butterflies similar to the case of the honeycomb lattice. However, 
in contradistinction to the honeycomb lattice, one notices the presence of a dispersionless 
band at E — 0, which is flat with respect to the variation of the magnetic flux, and is 
protected by a gap opened at B ^ (2J 13] . The spectrum exhibits Bloch-Landau bands at 
the extremities and also relativistic Dirac-Landau bands towards the middle. The two types 
of bands are distinguished by opposite chirality dE/dcj) and by different dependence on the 
magnetic field. 

The periodic boundary conditions discussed above can be properly used for describing 
infinite lattices, however when interested in mesoscopic plaquettes they have to be replaced 
with vanishing boundary conditions. We intend to identify the differences introduced by the 
finite size, which will turn out to be non-trivial in the case of the Lieb lattice. 

For the confined Lieb lattice, the eigenfunctions can be obtained as combinations of 
functions Eq.(9) or Eq.(10) with coefficients that ensure the vanishing of the eigenfunction 
along the edges. As a technical detail we mention that (along the Ox direction, for instance) 
the finite plaquette begins with the atom A in the first cell, and also ends with an atom A 
which belong to the (iV+l)-th cell. This means that the wave function \&(k) > should vanish 
at the site B in the 0-th and (N + l)-th cell, i.e.: < ®(k)\tf N+ ^ m \0 >=< $(£)|6jj0 >= 0. 
Similarly, the vanishing condition along Oy occurs at the site C in the 0-th and (M + l)-th 
cell along this direction, i.e.: < 3>(&)|4, M+1 \0 >=< $(^)|c^ |0 >= 0. 

In the localized representation, which is the proper one in the case of confined systems, 
the eigenfunctions > corresponding to E = Qq = look as follows: 

- JV+l Af+l r, , fc„ 



2 ™^ \ 2t v cos^ I 

T" 

2t x cos ! f . I 



$°(/c) >= \ — — t\ — — - ( , = 2 ^=sink x n sinkJm - -) 6* 10 > 

KJ V N + l V Af + l ^ ^ A / A| 2 + A 2 2 J nml 

n=l m=l VI I 1 I I 



, zsinkJn ) sink v m c^lO > ), (11) 

7|A| 2 + |A| 2 2 ; v nml h y ! 

where k x , k y are obtained from the condition that the wave function vanishes at the boundary, 
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and equal k x = pir/ (N + 1) (p=l,..,N+l) and k y = -^j-j- (q=l,..,M+l). Since the situations 
p = N + 1 (at any q) and q = M + 1 (at any p), generate |$° >= 0, we are left in Eq.(ll) 
with only NM non-vanishing degenerate orthogonal eigenfunction. 

The eigenfunctions |$ ± (/c) > corresponding to the other two energy branches can be 



written similarly as: 



'N+l M+l 



-> / 2 ' ^ v— -\ 1 1 
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2t x cos\ 1 . 
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H =sink x (n )sink v m c ' > , (12) 

v/|A| 2 + |A| 2 2 ; y nml ; ' v ; 

where states of the type A are this time also present. One can readily see that the number of 
non- vanishing states in each spectral branch is (iV + 1)(M + 1) — 1, since the point V = (tt, tt) 
has to be treated separately This is because its corresponding energy vanishes and the state 
should be counted in the flat band. In this case the wave function becomes: 

/ ; / ; N+l M+l 

|S° a >=: >= V^V^ EE!- 1 )"" 1 oL|0 > • (13) 

n=l m=l 

For the finite Lieb plaquette with vanishing boundary conditions, one may conclude that the 
flat band degeneracy equals NM + 1, while each other branch contains NM + N + M states, 
so that the total number of states equals indeed the number of sites 3NM + 2(N + M) + 1. 

In the presence of the magnetic field, the vanishing boundary conditions give rise to edge 
states which fill the gaps of the Hofstadter spectrum corresponding to the periodic system. 
Besides the edge states existing in the energy range of the Bloch -Landau levels (which are 
the only met for the finite plaquette with simple square structure), there are edge states 
in the relativistic range which show opposite chirality |16j . but also non-conventional edge 
states lying in the central gap which protects the zero energy dispersionless band. This 
last new class of edge states exhibits oscillating chirality when changing either the magnetic 
flux or the Fermi energy. These states will be studied in the next chapter. The fate of the 
zero-energy states in the presence of confinement will be discussed in the next section. 

The Lieb lattice can be generated from the simple square lattice by extracting each the 
second atom when moving along both Ox and Oy direction. Formally, this means either 
to push to infinity the energy E d of these atoms or to cut down the hopping integrals t' 



connecting them to the neighboring atoms, and it is instructive to follow the change of the 
spectrum when E^/t' — > oo. By driving the system in this way from 1 to 3 atoms/unit cell, 
the lattice periodicity is doubled along both directions, and the flat band is generated. The 
middle panel of Fig. 3 shows how the butterfly wings break off during the process giving rise 
to the relativistic range. 




FIG. 3: The energy spectrum as function of the magnetic flux for three values of the hopping 
integral t' (see text) : (left) t' = 1, corresponding to the simple square lattice, (middle) t' = 0.5, 
(right) t' = 0, corresponding to the square lattice with centered lines (Lieb lattice). <f> is the 
magnetic flux through the unit cell of the simple lattice measured in quantum flux units. The 
Hofstadter butterfly is obvious for t' = 1, while a doubled butterfly results for t' = in each of the 
intervals <j) £ [0,0.25], (j) 6 [0.25,0.5], etc. (one has to keep in mind that the flux through the Lieb 
unit cell is four times larger than (p). The energy is measured in units of hopping integral t. 

ILL SPECIFIC ASPECTS OF THE FINITE LIEB PLAQUETTE IN MAGNETIC 
FIELD: ZERO ENERGY FLAT BAND AND TWISTED EDGE STATES 

A. The properties of the flat band 

There are several pertinent questions which can be asked concerning the flat band in the 
energy spectrum of the Lieb finite system: what is the degeneracy, what is the response to 
the magnetic and electric field and to the disorder? 

Let us find first the conditions which should be satisfied by the zero energy eigenfunction 
\& e=o- Let H be the tight-binding Hamiltonian of a finite system and \1/ e its eigenfunctions: 

H = ^E n \n >< n\ + ^t nm \n >< m\, ^ = ^a n |n>, (14) 

n n,m n 



10 



where {\n >} is a basis of functions localized at the sites n. The condition e = 
generates a set of equations for the coefficients {a n }: 

E n a n + ^ t nraOL m = 0, Wl. (15) 
m 

Eqs.(15) are the necessary and sufficient conditions which must be fulfilled by the wave 
function in order to correspond to the zero eigenvalue E — 0. With E n = 0, and 
taking into account only the nearest neighbors (t nm = t) the above equations become simply 
Smev °- m = 0' f° r an y wnere t ne sum is taken over all sites in the first vicinity V n of the 
site n. In addition, if ^e=o anc ^ ^£=o are two degenerate states , the orthogonality condition 
reads J2 n a n a n = 0- The number of configurations {a n } which satisfy simultaneously the 
two conditions equals the dimension of an orthogonal basis in the space of the degenerate 
eigenf unctions at E = 0. 
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FIG. 4: (Color online) The three eigenstates of the flat band for a Lieb lattice composed of two 
cells. The eigenfunctions are ^^ ^ = ^2 nm ct nm \nm > and the coefficients a nm are indicated. We 
notice that the condition for the flat band appearance Snmev nomo a nm = holds for any site 
{n m }. 

An instructive illustration is the Lieb plaquette consisting of two cells (see Fig.4). The 
plaquette contains 13 atoms (6 of type A, 4 of type B and 3 of type C). There are three 
configurations of the coefficients a n which satisfy the conditions discussed above and they 
are pictured as (a), (b) and (c). (The numbers {0, —1, 1, —2} mentioned in Fig.4 represent 
the values, up to the normalization factor, of the coefficients a n ). 

With the notations used in the Hamiltonian (1), the three states can be written as: 
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0,q 


i = 0) = 
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+ 4i 


- CtL + a 12 - 4l + a 32] 1° > 
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0,q 


) = 0) = 
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4 - - 4 + b\ 2 + b\ 2 \ 10 > 


(16) 
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It is obvious that ^ n a l n = for any % — 1, 2, 3 and that < ^il^j >= for any i, j = 1, 2, 3, 
i.e. the three states correspond to E — and are mutually orthogonal. 

Next, we want to find out how the zero energy states Eq.(16) respond to a perpendicular 
magnetic field. In order to answer this question, we write the Hamiltonian (1) as: 

H{<t>) = # (O) (0 = O) + # (1) (0), 

= E i^mbnm + &L«»+l,m) " 1) + ff.C., (17) 

rim 

and perform degenerate perturbation with respect to ifW. Applying this approach to the 
two-cell Lieb system the matrix elements involved are < ty^lH 1 ^^ >= 8i siring and 
< ^^{H 1 ^^ >= and the secular equation reads: 



t -E 8i sinnd) ^ 



det 



0, 



-8i siniKp —E 
y -E j 

giving rise to the eigenvalues: E 1>2 = ±8tsinn(p and E 3 = 0. 

One remarks that the bulk state \J/ 3 does not couple to the magnetic field and its eigenen- 
ergy remains E 3 = 0. On the other hand, the surface states \J/i )2 get a dispersion which 
depends on 0. The conclusion of the perturbative calculation is that the magnetic field 
reduces by 2 the degeneracy of the zero energy band. 

Let us generalize now to a finite Lieb lattice containing N cells along the Ox-axis and 
M cells along Oy-axis, so that the total number of cells is N ce u = NM and the number 
of states is 3NM+2(N+M) + 1. It has been proved in the previous chapter that, at zero 
magnetic field, the number of zero energy degenerate states is N ceU + 1. Then, the two-cell 
model shows that in the presence of the magnetic field two states separates from the bunch 
so that the degeneracy of the flat band becomes N ceU — 1. Using a similar approach for 
the general case, one has to use the eigenfunctions Eq.(ll) and Eq.(13) and the expression 
Eq.(17) as the perturbation. One finds out easily that < $°(A;)|i/ (1) <I> (/?) >= 0, and that 
the only nonvanishing matrix elements are X{k) =:< $ (/c)|if >. In the general case, 
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the secular equation becomes: 

/ —E 




det 



-E 








\X(k x ) X(k 2 ) x 3 



... X{k 1 )\ 
... X(k 2 ) 

... —E ) 



which in the polynomial form reads E N ~ 2 (E 2 — X 2 ) = 0, where X 2 = X 2 (ki) + ...X 2 (k N ^i). 
This formula (where N stands here for the degeneracy of the flat band) says that from the 
whole bunch only two levels get a dispersion depending on 0, meaning that the degeneracy 
of the zero energy level is reduced by 2 in the presence of the magnetic field. So, the general 
finite Lieb plaquette behaves similarly to the two-cell model. 

The numerically calculated energy spectrum of the finite plaquette in perpendicular mag- 
netic field is shown in Fig. 5a, where one can check again the presence of the two levels which 
separates from the flat band while the most of the bunch at E — consisting of N ce u — 1 
states remains dispersionless. 




FIG. 5: The Hofstadter-type spectrum of a finite Lieb lattice of dimension iV? e 
for the clean plaquette, and (b) for the disordered one (disorder strength W 
measured in quantum flux units. 

The strong degeneracy of the flat band can be however lifted by a disordered potential. 
The broadening of the band depends on the strength of the disorder, however it remains 
independent of the magnetic field as in the case of the clean system (see Fig. 5b). We use a 
diagonal disorder of the Anderson-type characterized by the width parameter W [T7j. The 
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calculation of both the inverse participation number (IPN) and of the interlevel distribution 
indicate that in the middle of the disordered band the states are still delocalized, and 
described by the orthogonal Wigner-Dyson distribution ((3 = 1) which is the typical result 
in the absence of the magnetic field. This proves once more the absence of response of the 
flat band to the perpendicularly applied magnetic field, even in the presence of disorder. 
The inverse participation number (IPN) is defined as: 

IPN E = J2\< n \^E>\ i (18) 

n 

and indicates the degree of localization of the states. The small values of the IPN for 
energies in the middle of the density of states denotes the presence of extended states, and, 
as expected, the localization increases towards the band edges. The numerically calculated 
density of states and the dependence on energy of the inverse participation number are shown 
in Fig. 6a. Further information about the localization and the response to the magnetic field 
is provided by the distribution function of the level spacing between consecutive eigenvalues 
s n = E n — E n _i of the disordered system. Let us define the dimensionless quantity t n = 
s n J < s n >, where < s n > is the mean level spacing. In the disordered system, in the range of 
delocalized states, the level spacings are distributed according to the Wigner-Dyson surmise 

V {t) =h ^ e -^t\ (19) 

where /3 = 2 in the presence of the magnetic field, and /3 = 1 if I? = 0. As a signature of the 
distribution, the variance of the level spacing St =< Ss > / < s > is < St >= 0.4220 in the 
first case, and < St >= 0.5227 in the second one. Fig.6b shows the numerically calculated 
variance of the level spacing distribution, and one can notice that, in the middle of the flat 
band, where the states are delocalized, the variance is < St >= 0.5227. This means that, 
despite the presence of the magnetic field, the flat band behaves according to the orthogonal 
((3 = 1) Wigner-Dyson distribution instead of the unitary one ((3 = 2), as it is expected at 
B ^ 0. 

Another way to lift the degeneracy of the zero-energy band is to apply an in-plane static 
electric field. We expect specific aspects coming from the existence of the edges and of 
the lattice structure. In the numerical calculation the electric field applied along Oy-axis 
is simulated by replacing the atomic energies E nm with E nm + £y n , where y n is the site 
coordinate along Oy. Fig. 7a shows how the eigenvalues stemming from the flat band are 
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FIG. 6: (Color online) (a) The density of states and IPN in the flat band range as function of the 
energy for a disordered Lieb plaquette of dimension 20 x 20 cells averaged over 1000 configurations 
(disorder strength W = 0.3). (b) The variance of the level spacing distribution as function of 
energy for the same disordered system; the horizontal lines correspond to 0.4220 (as for the unitary 
ensemble) and 0.5227 (as for the orthogonal ensemble). 

split in several degenerate mini-bands which develop a Stark fan with increasing electric field. 
It can be checked that the number of mini-bands equals the number of lattice cells along the 
direction of the electric field. A perpendicular magnetic field gives rise to supplementary 
fine splitting and to the presence of states between mini-bands. This can be seen in Fig. 7b 
and also in Fig. 8. We have noticed that the flat band states are much more sensitive to the 
electric field than the edge states, and they give rise to a Wannier-Stark ladder at values of 
the electric field £ for which the edge states are still non affected. We have also numerically 
observed that the wave function in the I — th miniband is mainly localized in the I — th row 
of cells in the direction of the electric field. 

We already have seen that the flat band states do not show any diamagnetic response, 
and it is somehow surprising that the Wannier-Stark states coming from the former flat band 
exhibit a diamgnetic moment when the magnetic field is applied. It is interesting that each 
mini-band shows both positive and negative magnetic moments, and Fig. 8a suggests that 
the chirality dE/d(j) changes the sign at the center of the mini-band. We have studied also 
the localization properties of the eigenstates, particularly the localization along the edges 
P^ d9e , defined as: 

p E dge= J2 |$ a (*)| 2 , (20) 

i^Edge 
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FIG. 7: The low energy spectrum of a finite Lieb plaquette as function of the electric potential 
applied on the plaquette in the Oy-direction at: (a) 4> = 0, and (b) (f> = 0.12. 
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FIG. 8: The behavior of the flat band in crossed magnetic and electric field, (a) The orbital 
magnetization M a , and (b) the edge localization P^f 9e (b) vs. energy E a for a finite Lieb lattice 
of dimension N* ell = N^ ell = 15. The flat band turns into a set of 15 minibands, every miniband 
being composed of two parts with opposite magnetization. The states in the lowest and highest 
miniband have significantly increased edge localization (<p = 0.12 and e£L y = 0.2). 

where the index a indicates the state, and the sum is taken over all the sites which belong 
to the plaquette boundary. It turns out that the states which belong to the mini-bands from 
the extremities of the fan spectrum are strongly localized along the edges (Fig. 8b). The 
localization is of electric origin since the picture is similar no matter whether the magnetic 
field is present or not. 

We conclude, saying that the disorder lifts the degeneracy of the flat band keeping the 
states independent of the magnetic field, while the electric field produces states which re- 




16 



spond to the magnetic field and show specific diamagnetic moments. 
B. The twisted and type II edge states and their properties 

The confinement of the Lieb lattice induces several types of edge states. Besides the 
conventional edge states found in the Bloch-Landau and Dirac-Landau regions, there are 
still two other classes of edge states. We discuss first the twisted edge states lying in the 
magnetic gap opened around the degenerated energy level E = 0. Although the new states 
are localized along the perimeter of the plaquette, they do not follow the known behavior 
of the conventional edge states. The new class of edge states manifest specific properties: i) 
their energy depends on the flux in a periodic way. This means that the chirality defined by 
the sign of dE/dcf) is not conserved but alternate when changing the flux, in contradistinction 
to the usual edge states either in the Bloch-Landau or Dirac-Landau domain. Obviously, the 
alternate chirality should be reflected also in oscillations of the orbital magnetization at the 
variation of the magnetic flux, ii) their energies as function of the flux appear as twisted into 
bunches; for the clean square plaquette shown in Fig. 9a each bunch consists of four states, 
iii) the states prove the lack of robustness against disorder and iv) prove specific transport 
properties, namely, the twisted edge states carry a finite longitudinal resistance accompanied 
by vanishing Hall resistance. A piece of the spectrum of the clean plaquette in the energy 
range of twisted edge states is shown in Fig. 9a, where bunches consisting of four twisted 
edge states can be observed. One also has to notice that, at a given flux, the states in the 
bunch may show opposite chirality meaning that they carry diamagnetic currents moving 
in opposite directions. In the presence of disorder (Fig.9b) one notices that the twisted 
eigenenergies get stretched but the rest of the spectrum (the band and the edge states in 
the Dirac region) is not affected. This indicates that the twisted states are very sensitive 
to disorder. The understanding of this effect is simple in the sense that the degeneracy at 
crossing points [19j is lifted by the perturbation introduced by the impurity potential, and 
this occurs even at weak disorder. The Lieb lattice exhibits still another specific edge states 
(which we call type-II edge states), which in Fig. 12 are placed immediately above the Dirac- 
Landau bands at the transition from Dirac bulk to conventional edge states. They cannot be 
identified according to the sign of the magnetic moment |20j since their chirality dE n /d<f) is 
the same as for the bulk (band) states [IE]. Nevertheless, the diamagnetic currents of these 
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FIG. 9: The eigenenergy in the range of the twisted edge states vs. the magnetic flux eft for a 
pure Lieb lattice (a) and for a disorder Lieb lattice (b). The twisted edge states has an oscillatory 
behavior when the magnetic flux is varied, and they form bunches with four states in each bunch. 
The oscillatory behavior is destroyed by disorder in the right figure, but the conventional edge 
states (shown in the lower part of the spectrum) remain robust against disorder. The dimension 
of the Lieb lattice is N£ lt = ' ,, = 10 and the amplitude of the Anderson disorder is W=l. 



— 'S^— — ) • 

(a) (b) 

FIG. 10: The behavior of the edge states with disorder, (a) The absence of the backscattering for 
a conventional edge states. A pair of twisted states which are sufficiently close in energy may suffer 
the backscattering suggested in (b), which induces the localization shown in Fig. 11. 

states are located along the edges of the plaquette. These edge states show a double-ridge 
profile and carry current in both directions, but nonetheless the total magnetization remains 
of bulk-type. 

In Fig. 13 the diamagnetic currents of bulk states, type-II edge states and of conventional 
edge states are sketched. The twisted edge states may show currents similar to both con- 
ventional and type-II edge states. Compared to the twisted states, the type-II edge states 
behave substantially different in the electronic transport. These states will be studied in de- 
tail elsewhere. The contribution to the magnetization of each eigenstates \a > is calculated 
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FIG. 11: l^j 2 calculated for a conventional edge state in the Dirac range (left) and for a twisted 
edge state (right) for a disordered plaquette with W = 0.2. One observes that this low disorder 
does not affect the conventional edge state but localizes the twisted state. 




magnetic flux 

FIG. 12: The eigenenergies vs. the magnetic flux <f> for a pure Lieb lattice in the range which 
emphasizes the type II edge states. In the spectrum, they are located between the bulk states in 
the Dirac-Landau range and the conventional edge states of the first gap. Their energy decreases 
with the magnetic field similar to the bulk states, however they have edge localization of the wave 
function. The dimension of the Lieb lattice is TV? ' = ,, = 10. 
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FIG. 13: The sketch of the diamagnetic currents in the Dirac-Landau range of the spectrum: (a) 
the counterclockwise loop of a bulk state, (b) the double ridge current of a type-II edge states, and 
(c) the clockwise loop of a conventional edge states. The twisted edge states may show both (b) 
and (c) aspect. 



following the approach from [2T] , namely: 

M a = 



.dH. 
< ot — — a > 



int. 



inn 

-imw 



me 



< a\A nm >< B nm \a > -me™* < a\A nm >< B n _ hm \a > 



+ me- im ^<a\B nm ><A n+l ^\a>-me im ^<a\B nm ><A nm \a>). (21) 

All the matrix elements in the above equation are known once the eigenstates \a > are 
calculated numerically in the presence of the magnetic flux. Fig. 14 depicts the diamagnetic 
moments and the localization at the edges of different type of states. The bulk (band) states 
show positive magnetization and vanishing localization at the edge, the conventional edge 
states show negative magnetization, and 97% localization at the edge. The twisted edge 
states show alternating magnetization, as expected, but also an unanticipated differences 
in the degree of edge localization. This is because they exhibit either a single- or double- 
ridge wave function. Obviously, the double-ridge wave function is not strictly stuck to the 
edge so that P e d ge ~ 0.7 — 0.8, while for the single-ridge states the same parameter goes up 
to 0.9. Fig. 14 points out that the single-ridge states which are localized close to the edge 
exhibits negative magnetization (as the conventional edge states), while the double- ridge 
states exhibits positive magnetization. 

IV. THE INTEGER QUANTUM HALL EFFECT 



The quantum transport of the 2D Lieb plaquette shows some similarities with the case of 
graphene, however it also reveals particular properties. The Hall resistance as function of the 
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FIG. 14: (Color online) Magnetization M a and localization at the edges P e dge corresponding to 
the eigenenergies E a : the bulk (band) states (green) show positive magnetization and vanishing 
localization at the edges; the conventional edge states (blue) show negative magnetization and 97% 
localization at the edges; the type-II and the twisted states show 60 — 80% localization at the edges. 
The data are for a clean Lieb plaquette of dimension N£ ell = = 10 and <f> = 0.16. 

Fermi energy at a given quantizing magnetic field was obtained in Ref.3 by calculating the 

Chern numbers, and has the general aspect which can be observed in Fig. 15 which we obtain 

in the Landauer-Biittiker formalism: starting from the bottom of the spectrum, Rh shows 

h/e 2 steps in the Bloch-Landau region, then change the sign, and shows again h/e 2 steps in 

the Dirac-Landau region. The steps of the quantum Hall plateaus differs from those of the 

honeycomb lattice since in the Lieb case there is only one Dirac cone per the unit cell. The 

l 

change of sign is associated with the opposite chirality of the edge states in the two regions 
and occurs at E — ±2t, while in graphene the same change takes place at E — zkt [22]. The 
density of states (shown in blue in Fig. 15) is calculated as DOS = — - J2 n ^ m ^nni.^)i where 
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G + is the retarded Green function for the mesoscopic plaquette connected to the leads. In 
order to calculate the transport properties the mesoscopic plaquette must be connected to 
leads, the whole system being described in the tight-binding approach by the Hamiltonian: 

H = H s + H L + tH ls , (22) 

where the first term is just the Hamiltonian (1), the second term describes all the leads, and 
H LS couples the leads to the plaquette with the strength r. With G+p(E) =< a\(E — H + 
i0) _1 |/3 >, the electron transmittance between the leads a and (3, in the Landauer-Buttiker 
formalism, is given by: 

T a ,p = Ar'\Gl fj (E F )\ 2 Img^E F )Img L p(E F ), (23) 

where g L is the Green function of the isolated leads. In what follows we shall discuss the 
interesting question of the contribution to transport of the twisted edge states introduced in 
the previous section. The answer can be found in Fig. 15 in the energy range E G [—0.8, —0.6], 
where one observes that the twisted edge states found in that range do not support the Hall 
resistance (Rh = 0), however they contribute to the longitudinal resistance, which exhibits 
an oscillating behavior. It is also to notice that all the oscillations minima equals Rl = 1/4, 
a fact that should find its explanation. 

In exploring these unexpected effects, the first step should be to identify the transmittance 
matrix. The numerical investigation presented in Fig. 17a shows that, in the range of the 
twisted edge states the properties of transmittances T a> p are very specific: they are not 
quantized, show an oscillating dependence on the energy, and, mainly, satisfy the symmetry 
relation: 

T a ,a+1 — T a+ i t0l , (24) 

while in the range of the conventional edge states, where the quantized plateaus occur, the 
usual properties of quantum Hall effect hold: T a Q+1 =integer and T a+lja = (for a given 
sign of the magnetic flux). Combining Eq.(24) with the general property J2 a T a> p = 0, which 
expresses the current conservation, it turns out that the transmittance matrix for the edge 
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FIG. 15: (Color online) The transport properties of the Lieb lattice under the magnetic field: 
Hall resistance Rh, longitudinal resistance Rl and density of states DOS for a finite Lieb lattice 
connected to four transport leads. The quantum Hall effect can be observed for E G [—2.75, —0.8] 
(Rh integer and Rl = 0). In the Bloch-Landau part of the spectrum (E < —2) one has Rh < 0, 
while in the Dirac-Landau part (E > —2) one has Rh > 0. In the energy range E S [—0.8, —0.6] 
the transport properties are due to the twisted edge states and we get zero Hall resistance Rh = 
and oscillations of the longitudinal resistance with the characteristic minima at Rl = 1/4. The 
density of states exhibits maxima at the transition between the Hall plateaus and for the energy 
values where the twisted edge states appear. The dimension of the plaquette is 10 x 30 unit cells, 
the magnetic flux is <fi = 0.12. The resistance is in units h/e 2 , DOS in arbitrary units, and the 
energy in units t. 



transport in the domain of twisted edge states can be written as: 

{ -2T T T \ 
T -2T T 

T = 

T -2T T 
\ T T -2T ) 

The transmittance T a ^ relates the current through the lead a to potentials at the contact 
sites (5 as: 

i a = jJ2 T ^ ( 25 ) 



and, for the four-lead device considered in Fig. 16, in the Landauer-Biittiker formalism, the 
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FIG. 16: (Color online) The four-lead Hall device: illustration of the edge currents carried by the 
twisted edge states, for which the symmetry relation T Qja+ i = T a+ i tQ holds. 

Hall and longitudinal resistance are given (in units h/e 2 ) by: 

Rl = -Rl4,23 — (T24T3I — T2lT 3 4)/D, 

Rh = (-Ri3,24 — ^24,13) /2 = (T23T41 — T21T43 — T 32 Ti4 + Ti 2 T 34 )/2D, (26) 

where D = —AT 3 is a subdeterminant of the matrix T. By the use of the above equations 
and of transmittance matrix T, valid in the range of twisted edge states, one obtains im- 
mediately, a vanishing Hall resistance (Rh = 0) and the longitudinal resistance Rl = 1/4T 
(in units h/e 2 ). The minima of Rl observed in Fig. 17b correspond to the maxima of the 
transmittance, and, obviously, the value Rl = 0.25/i/e 2 expresses a perfect conducting one 
channel transport with T = 1. It turns out that, although carried by edge states, the current 
shows a dissipative character. The oscillations of Tyi and T21 in the domain of the twisted 
edge states (E e [—0.95, —0.75]) follow the similar oscillations of the density of states, while 
in the quantum Hall regime (E e [-1.5,-0.95]) the DOS is flat. 

Another interesting problem is the transition between the first and second plateau of the 
Dirac-Landau region (J5 € [—1.55, —1.5]) which is much wider than similar transitions in the 
Bloch-Landau region. The transition get a width which is due to the presence of the type-II 
edge states (observed in the spectrum Fig. 12 above the Dirac band), and is accompanied by 
oscillations of the transmittance (see T12 and T21 in Fig. 17a). We have to stress that in this 
energy range T12 and T21 are no more equal, and the numerical calculation suggests that 
T 12 = T21 + 1. It means that the symmetry Eq.(24) remains specific to the twisted edge 
states. 
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FIG. 17: (Color online) (a) The transmittances Tyi and Tn showing the quantized values Tyi = 1 
and T21 = in the range of the IQHE, non-quantized oscillating values Tyi = T21 in the range 
of twisted edge states E > —0.95, and Tyi = T21 + 1 in the range of type-II edge states E G 
[—1.55,-1.5]. (b) The Hall and longitudinal resistance: Rh (blue line) vanishes in the range of 
twisted edge states, while Rl (red line) exhibits oscillations with minima Rl = 0.25. The minima 
of the longitudinal resistance occur at the energies where the transmittance (black dashed line) get 
the maximum value Tyi = Ti\ = 1. The dimension of the plaquette is 10 x 30 unit cells, cj) = 0.16. 

In order to figure out a scenario for the vanishing Hall effect, we remind first that, in 
the range of the QH effect, all the (conventional) edge states responsible for the plateaus 
of the transverse magnetoresistance get a unique chirality determined by the direction of 
the magnetic flux. This results in a definite sense (say clock-wise) of the current such that 
T a , a +i— integer and T a +i )Q! = 0. The symmetry T a)Q+ i = T a+ i^ a which occurs in the range 
of the twisted edge states is characteristic to the absence of the magnetic field. So, this 
symmetry indicates a 'loss of influence' of the magnetic field followed by a vanishing Hall 
effect. As a support of this idea we note that the twisted edge states show alternating (clock 
and anti-clock) chiralities which allow for the transmittance in both directions, as sketched in 
Fig 16. This might be an heuristic explanation for the vanishing of the transverse resistance 
despite the presence of edge states. 

Let us discuss now the effect of the disorder. It is known that the conventional edge 
states are robust to disorder, whereas the bulk states (which form the Landau bands) are 
more sensitive, so that the IQHE of a disordered plaquette shows robust plateaus, and 
a broadened transition region between two consecutive plateaus. On the other hand, as 
we have shown, the disorder localizes easily the twisted edge states, changing in this way 
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FIG. 18: (Color online) The Hall resistance Rh (red) and transmittances Tyi (green) and T21 (blue) 
for a disordered Lieb lattice in the transition region from the first Hall plateau to the domain of 
twisted edge states. The figure shows that the symmetry relation T\i = T21 for twisted edge 
state transport holds also in the presence of disorder. The Hall resistance is Rh = 1 for the 
energies corresponding to the first gap with conventional edge states, and Rh = for the energies 
corresponding to the twisted edge states. The dimension of the plaquette and the magnetic flux 
are the same as in Fig. 17, the disorder amplitude is W = 0.5. 

their transmittance properties. Fig. 18 depicts the transmittances T 12 ,T 21 for a disordered 
Lieb plaquette compared to the same transmittances of the clean system. One observes the 
quantized values in the range of the conventional edge states followed, in the range of twisted 
states, by reduced values of the disordered transmittance which replace the peaks specific 
to the clean system. It is worth to say that the symmetry T a ^ a+ i = T a+ i tCt is preserved for 
each individual disordered sample, and, as a consequence, the Hall effect vanishes similar to 
the clean case. 

V. CONCLUSIONS 

We have found that the specific topology of the 2D Lieb lattice induces remarkable 
spectral and transport properties. Up to point there are similarities with the electronic 
energy spectrum of graphene in what concerns the presence of a Dirac-type cone at low 
energy, however, in addition, a macroscopically degenerated flat band occurs at the middle 
of the spectrum. The perpendicular magnetic field applied on a finite (mesoscopic) Lieb 
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plaquette opens a gap around the flat band, and we show the presence in this gap of a 
new class of edge states with alternating chirality (which we call twisted edge states). The 
flat band is insensitive to the magnetic field, however an in-plane electric field, and also the 
disorder, lifts the degeneracy. The electric field applied on a finite (mesoscopic) system gives 
rise to a Wannier-Stark fan composed of degenerate mini-bands, the number of them being 
equal to the number of cells along the direction of the field. 

The macroscopic degeneracy of the flat band is lifted by disorder, and the degree of 
localization and the level spacing distribution are studied. It turns out that not only the 
ordered flat band, but also the disordered one, does not feel the magnetic field; indeed, we 
prove that the level spacing distribution of the disordered system follows the orthogonal (/3 = 
1) Wigner-Dyson distribution, which usually describes disordered systems in the absence of 
the magnetic field. 

We calculate analytically the orthogonal eigenfunctions of the finite Lieb system cor- 
responding to the three spectral branches in the low energy range, both for periodic and 
vanishing boundary conditions. In this way we find also the degeneracy of the zero energy 
flat band which, in the periodic case, equals the total number of unit cells N ce u (except when 
both N* ell and N^ eU are even numbers, in which case the degeneracy increases to N ce u + 2), 
while in the case of the closed boundaries the degeneracy is N ce u + 1. A toy model composed 
of only two unit cells helps to understand the behavior in the presence of a perpendicular 
magnetic field. The perturbative calculation shows that 2 states of the flat band separate 
from the degenerated bunch, and belong to the class of twisted edge states. 

The eigenenergies of the twisted edge states depend in an oscillatory manner on the 
magnetic flux, i.e., show an alternating chirality, and contrary to the conventional edge 
states, the diamagnetic moment change the sign when the magnetic flux is varied. These 
type of edge states generated by the magnetic field are not protected by the broken time- 
reversal symmetry, and proves to get localized even at low disorder, when the conventional 
edge states remain robust. 

The transport properties are calculated by attaching leads to the finite Lieb system and 
using the Landauer-Biittiker formalism. The quantum Hall resistance looks similar to that of 
the graphene except the steps are equal to h/e 2 (instead of h/2e 2 ), however in the domain of 
the twisted edge states the properties become unconventional: the Hall resistance vanishes, 
while the longitudinal one shows oscillations which can be correlated with the oscillations 
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of the density of states (calculated in the presence of the leads). This behavior stems from 
the symmetry of the transmittance T a>a+1 = T a+1>a which occurs despite the presence of the 
quantizing magnetic field. The symmetry holds also in the presence of disorder. 
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